Quantum skyrmions and the destruction of long-range antiferromagnetic order in the 
high-T c superconductors La2 x Sr x Cu04 and YBa2Cu306+:r 
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We study the destruction of long range antiferromagnetic order in the high-T c superconductors 
La2- a: Sr a ;Cu04 and YBa2Cu306+;c. The CP -nonlinear sigma model formulation of the two- 
dimensional quantum Heisenberg antiferromagnet is used for describing the pure system. Dopants 
are introduced as independent fermions with an appropriate dispersion relation determined by the 
shape of the Fermi surface. Skyrmion topological defects are shown to be introduced by doping and 
their energy is used as an order parameter for the antiferromagnetic order. We obtain analytic ex- 
pressions for the skyrmion energy as a function of doping which allow us to plot, without adjustable 
parameters, the curves Tn(x c ) X x c and M(x) x x, for the two compounds, in good quantitative 
agreement with the experimental data. 



PACS number(s): 74.72.Bk, 74.25.Ha 

High-temperature superconductivity is by now well 
established to arise from doping quasi two-dimensional 
(quasi-2D) Mott-Hubbard antiferromagnetic insulators. 
The doping process initially produces the destruction 
of the antiferromagnetic ordering, giving place to a 
quantum spin-liquid disordered phase pj. The two 
best studied examples are I^-zSr^CuC^ (LSCO) and 
YBa 2 Cu 3 6+;!; (YBCO) for which the Neel ordered 
ground state at x — is replaced by a quantum dis- 
ordered state for x c w 0.02 and x c « 0.41 ||, respec- 
tively, at T = 0. 

It has long been recognized that the pure compounds 
are well described in terms of an S = 1/2 quantum 
Heisenberg antiferromagnet (QHAF) on a square lattice. 
The long wavelength spin fluctuations of the latter, on 
the other hand, are described by the 0(3) quantum non- 
linear sigma model (QNLerM) in two space plus one time 
dimensions |^,^) whose Lagrangian density is 
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where n is the order parameter field, subject to the con- 
straint n 2 = 1, and p s and c are respectively the spin- 
stiffness and spin-wave velocity. 

The nonlinear sigma model possesses classical topolog- 
ically nontrivial solutions called skyrmions whose energy 
is Ef = Att p s . For fully quantized skyrmions, on the 
other hand, there is a reduction of the skyrmion energy 
to half of the above classical value |q| 



E s = 27rp s 



(2) 



Having in mind that the skyrmion energy is reduced 
by quantum fluctuations, it is natural to expect it to 
be further reduced by the extra fluctuations introduced 
through doping. Furthermore, since in the framework of 



the NLctM the skyrmion energy is proportional to the 
ground state magnetization we can use it as an order pa- 
rameter for the quantum phase transition associated to 
the destruction of the antiferromagnetic state. 

In this paper we revisit the continuum model proposed 
in Q, to describe the doping process in YBCO at T = 0, 
and extend it for including also the case of LSCO. The 
model predicts the creation of skyrmion topological de- 
fects precisely at the dopant's positions, as has been pro- 
posed earlier || (see also S). The formation of such 
magnetic textures causes a reduction of the ground state 
magnetization. As a consequence, the skyrmion energy 
itself is lowered and eventually vanishes at the Neel quan- 
tum critical point. By computing quantum defect corre- 
lation functions at T — 0, we obtain analytical expres- 
sions for the skyrmion energy as a function of doping, 
E 3 (x), which allow us to plot the curves M(x) x x, for 
both LSCO and YBCO, in good quantitative agreement 
with experiment. Subsequently, introducing finite tem- 
perature and interlayer coupling, we obtain the antifer- 
romagnetic part of the phase diagram, namely the curve 
Tn(x c ) x x c - As we shall discuss below, our analysis 
is compatible with a picture in which the formation of 
stripes would occur in LSCO but not in YBCO. 

The doping process: The chemical modification of par- 
ent compounds of the high-T c materials here considered 
produces the introduction of holes in the Oxygen orbitals 
in the layered Cu02 planes. In what follows, we shall de- 
termine how the skyrmion energy is modified from (|^) 
due to the presence of such holes. We propose that the 
doped system can be described in terms of a bulk (T = 0) 
0(3) QNLcM coupled to four-component fermion fields 
with dispersion relation determined by the shape of the 
Fermi surface. These represent the holes doped into the 
in-plane O p-orbitals, whereas the nonlinear sigma 
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field represents the spin density of the active electrons 
of the Cu ++ ions. For describing the coupling of the 
dopants to the Cu ++ spins, it will be convenient to make 
use of the CP 1 language, where n = zja^zj, with oV, 
being the Pauli matrices and Zj, i = 1, 2, complex scalar 
fields. Then, following M (see also jlO|]), we minimally 
couple the fermions to the CP 1 fields. This is also con- 
sistent with previous results where a minimal coupling 
of fermions to CP 1 fields was obtained in the long wave- 
length regime of the spin-fermion model |TT| . Let us con- 
sider the cases of YBCO and LSCO separately. 

a) For YBCO we have an almost circular shape for 
the Fermi surface, which is centered at k = [55 , 
see Fig. [j]a. We can then use the dispersion relation 
e(k) = ^Jk 2 vp + (m*vp) 2 , with m* and vp being re- 
spectively the effective mass and Fermi velocity of the 
dopants. Observe that close to the Fermi level of the 
doped system [ep > m*vp) the above dispersion relation 
behaves as e(k) — ep ~ vp{k — kp) as expected. This dis- 
persion relation corresponds to a Dirac kinetic term for 
the fermions and the doped system can then be described 
by the partition function 

Z = J V[z,z,A^,^] 6[zz-l] ^"-A"] 

x exp j jf dr J d 2 x [2p a (D„Zi)* 



(3) 



where we allow for a residual Hopf term Ch = 
(9 /2)e Aia/5 ''ApdaAp, which has been shown to exist in the 
presence of topological defects on the ground state ]T3| . 
As usual, Dfj, = + iA^, with <9 M — (d T /c, V), and 
A^ = izjduZi is the Hubbard- Stratonovich CP 1 vector 
field. 





FIG. 1. Approximate Fermi surfaces for: a) YBa2Cu306+;c 
and b) Laa-^Sr^CuO^ f[§. 

The second delta functional constraint in (||) is used in 
order to introduce the in-plane hole concentration pa- 
rameter S. In its argument, j M = ip'y^ip and A M = 
4(5 f XL d^S 3 (z — £) for a dopant introduced at the posi- 
tion X and moving along the line L. The factor of four 
in the definition of A M accounts for the degeneracy of 



the adopted representation (4-component) for the Fermi 
fields. The zeroth component of the associated Lagrange 
multiplier will be the chemical potential. 

It has been shown in 0, that upon integration over 
the fields z,z,ip,?p, the resulting equation of motion for 
Ao is such that a skyrmion topological defect configura- 
tion concides with the dopant position at any time and 
ird = 26. In other words, holes dress with skyrmions 
and we see that indeed a Hopf term is required in (0) for 
nonzero doping. We stress that, because of the CP con- 
straint, a mass term is generated for the A^ field. There- 
fore, as a consequence of screening, there will be neither 
statistical transmutation nor the generation of spurious 
in-plane magnetic fields which would be inconsistent with 
muon spin relaxation experiments E3 . 

The parameter S counts the number of holes in the 
Cu02 planes. This must be connected to the oxygen 
stoichiometry parameter x. For YBCO, it is known that 
the out of plane O-Cu-0 chains play an important role 
in the process of doping. At low doping, x < 0.18, most 
of the holes go to the out of plane chains and the system 
can be considered as pure. In view of this we propose 
that the density of in-plane charge carriers is related to 
the oxygen stoichiometry by 8 = x — 0.18. 

For evaluating the skyrmion energy E s , we use the 
fact that skyrmions are topological defects whose quan- 
tum properties can be fully described by disorder field 
opertors p |l5p|]. Using these, we are able to calculate 
the large distance behavior of the skyrmion correlation 
function, ((J(X)n(Y)) -► e -( s -./ fic )l x - y l/|X-Y| I/ , from 
which we can extract the skyrmion energy E s . For the 
pure system we have that E s is given by (Q). For the par- 
tition function (ph, we obtain, after integration over z, z, 
and ip, %/} fields, Wf 



E S (S) = 2np s 1 
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where 7 = 32 ^ 9 ^ 16 ^ 6 ' > = 10.9398 is a numerical factor- 
that comes from the integration over the fermions, and 
api = a/s/2 = 2.68 A, the minimal distance between two 
oxygen atoms, is the lattice spacing for dopants. The 
above skyrmion energy can be put in the form E s = 
2-nps (S) if we define an effective 6 dependent spin-stiffness 
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We then conclude, after comparing (^) with (g), that 
the doped system can be described by a QNLerM with a 
generalized S dependent spin-siffhess given by (j^). 

We can immediately obtain the reduced sublattice 
magnetization as a function of doping as: M(x)/M(0) = 
v / p s (x)/p s . This is plotted in Fig. |for he = 1.00 ±0.05 
eV A and p s = 0.069 eV @. We see that our theoreti- 
cal prediction is in good agreement with experiment and 
p s (S) vanishes for 
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(6) 



Using the above experimental input, we obtain S c = 
0.23 ± 0.03 or x c = 0.41 ± 0.03 for the quantum criti- 
cal point. 
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magnetization for Yl^CuaOe+i. Experimental data from 
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b) For the case of LSCO, the Fermi surface is 
given approximately by the one of Fig. |l]b, which 
has also been observed in ARPES ]l7fl . We shall use 
the fact that the Fermi surface of LSCO can be ob- 
tained by shifting the one of YBCO towards Q = 
(■n / a,ir / a) and symmetry related points in the Brillouin 
zone. This leads to the new dispersion relation e(A:) = 
y/[(k x ± ir/a) 2 + [k y ± n/a) 2 ]v 2 F + (m*Vp) 2 , which cor- 
responds to a kinetic term for the fermions in (^|) modified 
by a shifted derivative term. Even though the partition 
function remains unchanged, this shift will have a non- 
trivial effect on the skyrmion correlation function and en- 
ergy because of the constraint on the fermionic density 
of states. The computation of the skyrmion correlation 
function proceeds as before [[7|(| and we now obtain the 
following expression for the skyrmion energy 



E s (5) = 2tt Ps 1- 
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There are two important differences between ([?]) and 
The first one is a linear S dependence in (^) which 
is absent in (Q). This was generated by the shift in the 
dispersion relation and is a consequence of the fact that 
the Fermi surface of LSCO is not centered at k = in 
the Brillouin zone. The second one is an extra factor of 
four multiplying 5 in (0), which is needed to account for 
the four branches of the Fermi surface of this compound. 

The above skyrmion energy suggests, in accordance to 
what has been done for the case of YBCO, that the doped 
system could presumably be described by a QNLerM with 
stiffness 
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However, the magnetization that would follow from this, 
according to the same procedure adopted for the case of 
YBCO, has a faster decrease with doping than the one 
observed experimentally, as we can infer from the solid 
line in from Fig. |^. This can be interpreted as a sign of 
stripes formation in LSCO Following we ob- 

serve that the presence of stripes produces an anisotropy 
which allows us to write p' s (S) = \J p%(5)p v s (5), where 
Ps(S) 7^ pl(S). Assuming that in LSCO stripes are par- 
allel to the y-axis, as the experimental results in |2(]] 
suggest, we propose that only the x-component of the 
spin-stiffness would be affected by doping, implying that 
p v s = p s and that p*(5) is given by (^). The resulting 
NLerM has an effective spin-stiffness 



p' s (S)=p s Jl- 
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Now the corresponding sublattice magnetization does 
agree with the experimental data as can be seen from 
the dashed line in Fig. |[ Both plots in this figure are 
obtained by using the experimental input Tic — 0.75±0.03 
eV A and p s = 0.051 eV Q. We have also used the well 
known fact that for LSCO, x = 5. Interestingly, the loca- 
tion of the quantum critical point is not affected by the 
formation of stripes. In fact, both p s (8) and p' s (S) vanish 
at 



1 J (2^/2~Tic/a) 2 + 4-fhcp s /ira D - (2y/2Tic/a) 

5 °=a — s ■ ( 10 ) 

4 Z^nc/naD) 

which, for the above experimental input, yields the value 
x c = 8 C = 0.020 ± 0.003 for the quantum critical point. 
We also remark that using the effective spin-stiffness (||) , 
implied by the stripes picture, the reduced sublattice 
magnetization M(x)/M(0) = \J p' s (x)/p s can be written 
in the form M(x)/M(0) = (1 - x/x c ) 1/4 , thereby pro- 
ducing a critical exponent of 0.25, which is very close to 
the value of 0.236, empirically obtained by Borsa et al. 
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The phase diagram: In order to obtain the critical line 
for the destruction of AF order in the Txx phase diagram 
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we must consider a nonzero interlayer coupling Jj_. For 
this purpose, we shall use the partition function pjj ] 



Z = J Vw.i exp | - 
+ (Vn,) 2 +«(n. 
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where a = (1/a 2 ) J±/J\\, with Jj_ and Jn being re- 
spectively the interlayer and intralayer couplings of the 
underlying microscopic quasi-2D QHAF and (3 = 1/T 
(ks = 1)- Contrary to the strictly 2D case now a finite 
value for the Neel temperature is obtained to order l/N 
in a large N expansion [El[ : 



T N = Airp s 
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Now, if we replace in ( 12j) the spin-stiffness p s by our 5 de- 
pendent expressions (|) and (||), respectively for YBCO 
and LSCO, and using J±/J\\ — 5 x 10 -5 , we obtain the 
phase diagrams plotted in Figs. |] and pj. 
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FIG. 4. Antiferromagnetic part of the phase diagram for 
YBa2Cu306+a;. For x c < 0.18, we assume a fixed Tjv = A2QK. 
Experimental data from |3|. 
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FIG. 5. Antiferromagnetic part of the phase diagram for 
L^-zSr^CuCU. Experimental data from |2|. 



In summary, our results indicate that the antiferro- 
magnetic phase of the high-T c cuprates can be correctly 
described in terms of a generalized NLcrM with an ef- 
fective, doping dependent, spin-stiffness which carries all 
the information about the quantum fluctuations intro- 
duced by the dopants. This can be inferred from the 
quantum skyrmion energy which is evaluated after inte- 
gration over the dopants in a model without adjustable 
parameters. Also, the results point towards a picture in 
which charged stripes presumably occur in LSCO, lead- 
ing to a slower decrease of the sublattice magnetization 
as a function of doping in comparison to the case where 
they are absent, see Fig. [|. Apparently, the only ef- 
fect of such phenomenon is to divide by two the quan- 
tum critical exponent of the sublattice magnetization (or 
the spin-stiffness) without affecting the location of the 
quantum critical point. Conversely, our results indicate 
that such phenomenon is absent in YBCO as our model 
(0) was able to correctly reproduce the zero temperature 
data for M (x)/M(0), see Fig. ||. This is consistent with 
density matrix renormalization group calculations in the 
framework of the n-leg ladder t — t' — J model |2|] . 

As a final comment we would like to remark that for 
the case of Bi2Sr2CaCu208+2; we expect that our result 
(^) shall correctly reproduce the experimental data with- 
out the scaling modifications done for LSCO. Indeed this 
material, while having a Fermi surface similar to LSCO 
p3[ , docs not seem to present the formation of stripes. 
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